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FORMULATION  OF  NORMALIZED  NON-LINEAR  FREE  ELECTRON 

LASER  EQUATIONS 

I.  Introduction 

Free  electron  lasers  ( FELs )  show  great  potential  for  becoming  a 
new  class  of  efficient  devices  capable  of  generating  intense  levels  of 
coherent  radiation,  continuously  tunable  from  sub-millimeter  to  beyond 
the  optical  regime. ^1-35^  The  FEL  is  characterized  by  a  wiggler  (pump) 
field,  for  example  a  spatially  periodic  magnetic  field,  which  scatters 
off  a  relativistic  electron  beam. 

In  this  paper  we  present  an  extension  of  the  formulation  of  the 
1-D  non-linear  self-consistent  FEL  equations  in  the  steady  state  amplifying 
configuration  with  space  charge  effect  for  a  circularly  polarized  wiggler 
field. 6 ,27 ’ 3 1 ^  We  included  all  the  efficiency  enhancement  schemes: 
i)  a  spatially  varying  wiggler  amplitude  and  period^26”28’31^  and  ii)  a 
D.C.  electric  field  EDC(z)  *  -  &zj^z  •  ( 3  2 ' 3  O  We  show 

that  all  the  efficiency  enhancement  schemes  are  somewhat  equivalent. 

The  complete  FEL  equations  cover  all  operating  regimes. 

For  a  variable  period  magnetic  wiggler  field,  the  FEL  equations  can 
be  appropriately  normalized  to  contain  only  three  parameters.  All  operating 
regimes  are  covered  by  the  three  parameters.  The  condition  for  space  charge 
to  be  neglected  is  derived.  In  the  absence  of  the  space  charge  effects, 
the  three  parameter  set  of  equations  is  reduced  to  a  two  parameter  set. 

In  the  high  gain  regime,  we  present  graphs  of  growth  rate,  efficiency 
and  saturation  amplitudes  of  the  radiation  field. 

In  the  low  gain  regime,  we  obtained  the  linear  gain  expression 
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with  the  space  charge  effect.^27’37*  38^  We  give  an  example  of  a  FEL  experi¬ 
ment  where  space  charge  effect  is  somewhat  important.  We  compared  results 
of  linear  theory  with  results  from  non-linear  calculations. 

II.  Normalized  FEL  Equations  with  Space  Charge  Effects 

We  present  a  self-consistent  non-linear  theory  of  the  FEL  in  the 
steady  state  amplifying  configuration,  see  Fig.  (1).  Included  in  our 
analysis  are  space  charge  effects  and  the  following  efficiency  enhance¬ 
ment  schemes:  i)  contouring,  spatially,  the  amplitude  and/or  the  wave¬ 
length  of  the  magnetic  wiggler  field,  and  ii)  by  applying  an  external 

D.C.  electric  field,  Enr  =  -  ^pc  ^ 

6z  z  ’ 

The  right  handed  circularly  polarized  external  magnetic  wiggler 

field  in  terms  of  the  vector  potential  is 

A0(z)  =  A0(z)(cos(  f  k0(z‘  )dz'  )&+  sin(  /*k0(z'  )dz*  )^  )  (1) 

*/o  *0  y 


where  kQ(z)  is  the  spatially  varying  wiggler  wavenumber.  The  vector  and 
scalar  potential  corresponding  to  the  radiation  and  space  charge  field  are 
represented  by 

A(z,t)  =  A(z) (cos(  f k+(z'  )dz'  -  u>t)^x-sin( Jk+{ z'  )dz'  -  cot}£M  ,  I 

*0  *0 
z 

(f>(z,t)  =  <Mz)cos( /0(k+(z'  )+k0(z'  ) )dz'  -  iot) 

+<j>2(z)sin(  f0(k+(z'  )+k0(z'  ))dz'  -  u>t)  ,  ( 

where  A(z),  <p  (z),  <j>  (z)  and  k,(z)  are  the  spatially  slowly  varying  amplitudes 
12  + 

and  wavenumber  of  the  fields.  The  frequency  ui  of  the  radiation  and  space 

charge  wave  is  taken  to  be  constant. 

Equations  governing  the  spatial  evolution  of  the  potentials  A(z,t) 

and<f(z,t)  have  been  derived  elsewhere. ^ 2 6"2 7 ^  We  begin  our  analysis  with 

the  fully  non-linear,  self-consistent  equations  for  A(z),  k+(z),  <^(z)  and 

4>  (z): 
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FEL  CONFIGURATION 


- INTERACTION  REGION - 

Fig.  1  —  Schematic  of  the  Free  Electron  Laser  model.  The  unmodulated  electron 
beam  enters  the  interaction  region  from  the  left.  In  this  figure  the  wiggler  (pump) 
field  builds  up  adiabatically  and  reaches  a  constant  amplitude  and  wavelength  for 
z  >  0.  The  pump  field  may  in  general  have  a  varying  period  and  amplitude,  which 
is  not  shown  in  this  figure. 
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> r*I- 


(w2/c2-  k2(z))A(z)  =  F  UiL  m0v  ^  F 
+  em  2c2  20 


A0(z)cos  <p(z,ta)  +  A (z); 


|a0 (z)sin  <Mz,t0)|  s 

2  2  V  „  /•  2^/  W 

(k+(z)  +  k0(z))4>z (z)  =  -  Fsc°~j2 — ^ J  dt0cos<Ji(z,  t0)  , 


^,2  2  V  f  2 Tr/w 

(k+(z)  +  k0(z))<n2(z)  =  -  F$c  ^r-—  dto  simp  (z,  t0)  ,  (3d) 

ip(z,  t0)  =  f  (k  (z1 )  +  k0(z‘  )  -  u>/v  (z1,  t0))dr  -  u>t0  ,  (4) 

*0  L 

where  wb  =  (4ir|e|  2n0/m0)^,  nois  the  local  beam  density,  Fgm  is  the  filling 
factor  for  the  electromagnetic  field,  Fsc  is  the  filling  factor  associated 
with  the  space  charge  field,  vZQis  the  initial  axial  beam  velocity,  vz(z,t0) 
is  the  axial  velocity  at  position  z  of  a  particle  which  was  at  z  =  0  at 
time  t0,  pz(z,t0)  =  Y  (z,to)movz(z,t0)  is  the  corresponding  axial  particle 
momentum  and 

y(z,t0)  =  (1  +  j e j 2 ( A0 ( z )  +  A(z,t  ))2/m?0c4+p2(z,t0)/m20c?)1'2  ,  (5) 

z 

where  T(z,t0)  =  /  vz(z|  tQ)dz'  +  tQ,  is  the  time  it  takes  a  particle  to 
reach  the  position  z  if  it  entered  the  interaction  region  z=0  at  time  tQ. 

The  phenomenological  filling  factors,  Fgm  and  Fsc  are  less  than  unity  and 
in  general  not  equal.  Roughly  speaking,  Fgm  is  the  ratio  of  the  cross 
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sectional  area  of  the  electron  beam  to  the  effective  area  of  the  radiation 

beam.  Since  the  radiation  field  is  not  directly  supported  by  the  electron 

beam,  Fgm  can  be  much  less  than  unity.  The  space  charge  field,  on  the 

other  hand,  propagates  on  the  electron  beam  and  falls  off  rapidly  away  from 

the  beam.  Hence  F  can  be  close  to  unity.  For  a  fixed  electron  beam  radius, 
sc 

F  decreases  as  the  radiation,  beam  radius  increases  while  F„  approaches 
em  sc 

unity  as  the  radius  of  the  cylindrical  waveguide  containing  the  electron 
beam  increases.  Physically,  the  reduced  beam  plasma  frequency  F$c  is 
less  than  the  local  plasma  frequency  u^.  The  parameter  F$c  can  also  be 
referred  to  as  the  plasma  frequency  reduction  factor.  In  the  usual  case, 
where  Fg[n  is  much  less  than  unity,  the  space  charge  wave  filling  factor  is 
close  to  unity,  FS(.al,  Chapter  8  of  Ref.  (31). 

Note  that  in  Eqs.  (3)  the  electron  beam  has  been  assumed  to  be 
initially  mono-energetic.  Generalization  to  an  initially  thermal  beam 
is  straightforward  and  considered  in  a  previous  one-dimensional  formulation 
of  this  problem. 

The  integrands  in  Eqs.  (3)  contain  the  particle  dynamics  which  are 
determined  from  the  Lorentz  force  equations.  It  is  straightforward  to 
show  that  the  axial  particle  velocity  v  (z,t0),  expressed  in  terms  of  the 
independent  variables  z  and  t0,  satisfies 


d  vz(z,t0) 
dz 


lei2  1  r  3  Ao  (z)-  3Ao(z)A(z)  ,  x 

ift1’  yU  3z  c°wU,t„) 


*A„(z)A(z)(k+(z)  ♦  k,(z)  -  _a-)sj,i  *(2,1.) 


+  Mz))(<Mz)  sin  <p( z,t0) 

rz'z,V  lel 

-*2(z)  cos  *<z.t.))]  *  iv^nr.TJc  ^  ' 


(6) 
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where  ya0=  (1  +  ( I ® I A0/m0c2 ) 2 J15,  yz(z,  t0)  =  (1  -  v|(z,  tj/c2)'*5  and 
Y(z,t0)  =  Yj.oYz(z,t0). 


From  (6)  we  can  derive  an  equation  for  the  phase  i p(z,  t0).  We  first 
define  vz(z,  t0)  =  vZQ  +  6vz(z,  t0),  Ak(z)  =  (k+(z)  +  k0(z)  -  u)/vz0)  and 
w  =  co0+Aw  where  |6v  |«v  ,u0  =  (1+8,J  y2  v  k„(0)  and  |Aw|«w0.  Noting 

L  l O  20  ZO 

that  <5 vz(z,  t0)  »  v2q  (3^(z,t0)/az  -  Ak(z) )/uj  and  k+(z)  +  k0(z)  -  vz(z,  t0)w/c2 
«2k0(0)(l  +  Ak(z)/k0 (0)  -  (2k0 (0) ) “ 1 a^(z ,  t0)/3z)  we  find  that  Eq.  (6)  takes 
the  form 

aVz,  to)  =  Mz)  2 1 e | 2  k0  (0)  d  A2  (z) 

az7  c4  ~  hi(z,  t0)  ~5z 


+  lelYo.L2  gj/c3 
moY03hMZ,  tfl) 


d4>cc 
~~dz  + 


4jeJi 

moc1* 


k20  (0) 

yIoH2  (z,  t0) 


A0(z)A(z) 

Yo 

h(z,to ) 


(1  1  ^  (z,tj 

h2(z,  to)  +  2k7Toy  az 

m  p  2 

je|-  (4>j (z)  sin  <P(z,  tQ)  -  <*>2 (z)  cos  i*/(z,  tQ)) 


^  sin  n,(z,  t0) 


,  (7) 


where 


h(z,t0)  =  yz(z,  t0)/YZo= 


^(z,t0) 

x  Ak(z) 

az 

ko  (0) 

In  obtaining  (7)  we  have  made  the  added  assumptions  that  A(z)  changes  little 
over  a  wiggler  period.  The  various  efficiency  enhancement  schemes  ((1)  con¬ 
touring  the  wiggler  period,  (2)  contouring  the  wiggler  amplitude  and,  (3) 
application  of  a  D.C.  electric  potential)  are  represented  by 
the  first  three  terms  of  Eq.  (7)  respectively.  The  affect  of  the  various 
schemes  on  the  phase  are  equivalent. 
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In  the  following,  we  will  consider  a  simplified  example 
where  the  D.C.  electric  field  is  zero,  =  ^DC  =  0,  the  am- 

DC  am¬ 
plitude  of  the  magnetic  wiggler  Aq  is  constant,  and  kQ(z)  = 

kQ(0)(l  +  e(z)).  We  will  show  that  the  FEL  equations  can  be 

normalized  in  such  a  way  as  to  be  functions  of  only  three 

parameters,  and  to  be  functions  of  only  two  parameters  in  the 

absence  of  space  charge-effects.  For  Y0>>1,  the  complete  set 

of  coupled  equations  in  (3)  and  (7)  can  be  written  as 


+2ttDS 

^3(z,  t0) 


<cos  ip  (z,t0)>sini|>  (z,  t0)  -<sin ^  (z,  t^)>cos<Hz 


»  t0)  > 


(6k(z)  -  e(z) )  a(z)  =  -  2n  <h"  (i,  l#)  costp  (z,  tQ)>» 

=  -  2Tr<h“l  (z,  t0)sin<J«{z,  t0)>  . 
dz 


(8a) 

(8b) 

(8c) 


where  z  -  kQ(0)z  is  the  normalized  axial  distance,  t0  =  u,t0,  6k(z) 

=  Ak(z')/k0(0)  +  Aw/uu0>  £k(z)  =  (k+(D  +  kQ(z")  -  w/v,  J  ,  a  =  u,0+Aw, 
=  (1+ezo)  Vz20vzoko(°^  I  Att'l<<“'o'  hl=,tn  )=  v7  f  z  ,  t  )  V-  0 

(l+Ak(z)/k#CO')-dit(z,ro)/ dz)  ^  .  <(...)>  =  f  (dt  0/2n)  (.  .  . ) 

^  2  rs^ 

is  the  ensemble  average  operator  a(T)  =  8ti(yz0  /  ij  A  (?)  /  A  0  *is 

the  normalized  radiation  field,  t,  =  \fiPm  i  V^7ck  (®))  is 

2  2  2 

the  beam  strength  parameter,  C0  =  (|e|A  5/ Vomoc  )  /2n  = 

7  2  2  2  2 

*tiS72n.  Do  =  (Fsc/FemH  / nYz 0  is  the  space  charge  parameter, 
c^ox  =  lelA0/Yorooc  *s  t*le  transverse  particle  velocity  due  to 


the  magnetic  pump  field. 


In  obtaining  Eq.  (5)  from  (3)  and  (4)  we  made  use  of  the 

2  2  2 

approximations  |  (A/2k+)  dk+/ dz  |  <<  dk/bz,  u-  -c  k+(zp  2uu 
0-ck+(z)),  vz(z,t0)  ^vzo+vzo  (d#(z,t0)/az  -  Lk(z))/ui  and 
where  appropriate  (k  (2)  +  k  (z))  »  w/c 

+  n 


The  space  charge  potential  is  given  by 


<Pi,z(z)  *  "  \  lip  p0-  Dc 
11  'zo 

where  Yx  0  =  y0/y,0=  (1  +  |e| 2k20/m20clt)l'2‘ 


<cosip  (z,  tQ)> 
<sin-4>  (z,  t0)> 


The  coupled  set  of  normalized  equations  in  (8)  completely  describe 
both  the  linear  and  non-linear  behavior  of  the  FEL  process  for  a  varying 
period  magnetic  pump  field.  These  equations  include  the  effects  of  space 
charge  in  both  the  low  gain  limit  (interference  regime)  as  well  as  the 
high  gain  limit  for  arbitrary  magnetic  pump  amplitudes.  Space  charge 
effect  is  manifested  in  the  last  term  of  the  modified  pendulum  equation, 
i.e.,  Eq.  (8a).  The  general  initial  conditions  needed  to  solve  the  FEL 
equations  in  (8)  are  simply  5k(0)  =  0,  e(0)  =  0,  h(0,t0)  =  1,^(0, t0) 

=  -  t0  and  3<j)(z,  t0)/9z|  =  -  Aw/w0  where  Aw  =w-w0  is  the  frequency 

mismatch  and  u>0  =  (1  +  B,  )yI  v,  k  (0). 

ZOZoZOO 

For  completeness,  the  linear  dispersion  relation^"^  can  also  be 
written  in  terns  of  the  normalized  variables. 

6k(<Sk  -  ?,)(6k  -C?)  *-wr*(l  +  6k)  (10) 
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In  the  absence  of  space  charge  effects  we  have  a  two  parameter  set 
of  equations  for  both  the  low  gain  and  high  gain  limits.  In  the  low  gain 
limit  the  two  input  parameters  are  C*  a(0)  and  Aw/o)0  while  in  the  high 
gain  limit  the  two  parameters  are  C*  and  Aw/co0.  All  operating  regimes 
of  the  FEL  are  therefore  completely  covered  by  two  parameters  when  space 
charge  effects  can  be  neglected.  The  inclusion  of  space  charge  effects 
introduces  a  third  parameter  D„.  The  importance  of  the  normalized 
equations  to  only  two  parameters  is  that  the  two  parameters  can  be  used 
to  identify  all  the  equivalent  FEL  experiments.  The  normalized  equations 
provide  a  better  understanding  of  the  trade  off  among  all  the  FEL  design 
variables. 


III.  Low  Gain  Regime  with  Space  Charge  Effects 

We  will  now  consider  the  low  gain  or  interference  regime  without 
applying  any  efficiency  enhancement  schemes.  We  may  expect  that  the 
single  pass  efficiency  will  be  small  which  implies  that  we  may  set  h  (z,t0) 
=  Yz(z»  t0)/yZ0  equal  to  unity  without  loss  of  accuracy.  The  fact  that 
h(z,  tjal  further  implies  that  both  [ Ak(z)/ko(0)  (  and  |3i|j(z,  t0)/3zj  are 
much  less  than  unity.  With  these  assumptions  the  non-linear  FEL  equations 
in  (8)  become 


4^1^-  +  C*a(0)  sin  ip  (2,  t0) 


+  2rr  Do(<cos^(z,  t0)>sinip  (z,  t0)  -  <sin  ^  (2,  t0)>  cos  ^  (z,  t0))  ,  (12a) 


6k(z)a(z)  =  -  2tt<cos^(z,  t0)>  , 

-jp  =  -  2Tr<sin^(z,  t0)> . 

Since  we  are  considering  the  low  gain  regime,  a(z)  is  replaced  by  a(0) 
on  the  right  hand  side  of  (8a). 

We  seek  for  solutions  where  space  charge  effect  is  moderately 
important.  To  solve  for  the  linear  dynamics  we  expand  i|»(z,  t0)  by  setting 
ip(z,  t  )  =  <Mz.  t0)  +  t0)  +  i^2^(z,  t0),  where  ii/^> 

The  equations  we  will  solve  are 


(12b) 

(12c) 


3Vo)(z,t0) 


=  0  , 


(13a) 


^  =  C2  a(0)sin  ^0^(z,  t0)  , 

9z 


3 ^ 2  ^ ( z »  t.) 


d<5k(2) 

dz 


(13b) 


(13c) 


+  2 nD*  (<cos  +  ^^)>sin  ^°^(2,  t0)  -  <sin(<J^°^+<f^  ^)>cos  i/°^(z,t0). 
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(14a) 


The  solutions  are 


t|/°)(z,t0)  =  t0  +  UZ  . 


^^(z,  t0)  =  -  Co  ^^-(sin(t0+y2)  -  sin  t0-  uz  cos  tQ)  , 

u2 

^  2  (z,  t0)  =  °  71  (2  sin  u z  -  uz  cos  uz  -  uz) 

u 


(14b) 

(14c) 


—  7T 


c2  On  a ( 0) 


U2Z2  ~  L13Z3  ~  ) 

sin(t0+  uz)  -  sin  t0-  uzcostG+  M  g  sin  t0+  — — cos  t0> 


where  u=  -  Auj/w0 

Substituting  (14)  into  (12c)  and  integrating  over  z  we  find  that 
the  gain  G  with  space  charge  effects  is  given  by 

G(e>0p)  =  ~7T4d|—  S(0,0p)  ,  (15) 

where 


s(e,ep)  =  2 |ep  ^  +  -jr  0p  ~w] 

9  =  yz/2  and  6  =-•  it2  D0z/2  /Fsc/P~m  £z/2yzo  is  the  contribution  due  t0 
space  charge  effects.  For  0p<  0,  we  can  rewrite  (15)  as 


S(0 


,0  )  =r  sin(0  +  6p)  "1 2  f  sin(0-0r)l2 

p  [  9 + 0p  J  '  L  6-0p  J  ' 


(16) 


When  space  charge  effects  are  negligible,  the  maximum  gain  occurs 
at  6=  1.3  where  a  (sina/0)2  /ae  =  0.54  and  is 


Wz>  -pT1-]  i  w) 

We  can  estimate  the  level  of  electromagnetic  field  needed  to  achieve 
saturation  in  an  interaction  length  L.  Electron  trapping  and  thus  satura¬ 
tion  occurs  when  the  phase  shift  of  the  particles  with  respect  to  the  wave 
is  approximately  it.  Loosely  speaking,  saturation  occurs  when |^^x^(2,t0)|55s 
7t ,  using  (12b)  we  find  that  this  condition  gives  a(0)«  2tt/ (C|k|L2 )  or 


A(0) 


(18) 
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•~T  (3 


o^zokoL>Z 


This  value  of  the  radiation  field,  necessary  to  achieve  saturation 
at  z  =  L  is  to  within  a  numerical  factor  of  order  unity,  the  same  as  that 
obtained  in  Ref.  (24)  using  a  different  line  of  reasoning. 

The  gain  expression  with  space  charge  effect  is  valid  when  l»e2, 
in  order  that  the  assumption  is  satisfied.  Therefore,  the 

gain  expression  is  applicable  when  the  space  charge  effect  is  only  moderately 
important.  The  condition  for  space  charge  to  be  important  is 


(  sM°)z\ 

v 2Y-  ) 


»  i 


(19) 


where  £  =  A'em  w0/  (/y^  ck0  (0) ) .  For  Fsc  =  Fgm  =  1,  this  condition  is 
identical  to  that  derived  in  Refs.  (37)  and  (38).  Condition  (19)  is 
dependent  on  the  length  of  the  interaction  region  L.  Space  charge  effect 
becomes  important  when  L  >(2yzo/£ko  (0) )  (F^/F^)*5 


The  more  appropriate  condition  for  neglecting  space  charge 
effect  is 


<  1 


(20) 


The  physical  interpretation  of  Eq.  (20)  is  that  space  charge  effects 
become  important  when  the  ideal  maximum  amplitude  of  the  space  charge 
potential  becomes  much  larger  than  the  initial  amplitude  of  the 
ponderomotive  potential.  This  condition,  rewritten  in  terms  of  the 
electron  density  is 

£  (koTz0A°A)  (2%m0c2Fsc)_1.  (21) 


1 
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The  condition  for  neglecting  space  charge  effects  Hq . 

(20)  is  not  directly  related  to  the  length  of  the  interaction  region, 
but  dependent  on  the  amplitude  of  the  vector  potential  of  the  radiation. 
For  the  low  gain  case,  using  Eq.  (18)  as  the  approximate  initial 
amplitude  of  the  radiation  vector  potential  for  obtaining  saturation 
at  length  L,  Eq.  (20)  reduces  to  Eq.  (19)  to  within  a  numerical  factor 
of  order  unity. 

The  linear  gain  expression  should  always  be  used  with  discretion. 
The  linear  gain  expression  is  valid  before  saturation  effects  becomes 
important.  When  the  frequency  mismatch  is  small, 


0<2^>  <  | 
y  o 

trap 


(22) 


some  particles  can  be  initially  trapped  in  the  ponderomotive  well, 
full  width  of  the  trapping  potential  is  given  by 

|e|*i 


The 


*L-\ 


trap 


vv 


-  2^a(0) 


(23) 


When  particles  are  initially  trapped,  saturation  effects  become 
important  roughly  within  a  distance  of  half  of  a  bounce  wavelength. 

A  bounce  wavelength,  L^,  is  the  distance  the  trapped  electrons 
travel  in  the  axial  direction  while  executing  a  bounce  period  in  the 
ponderomotive  potential  well.  Considering  only  nearly  resonant  particles 
in  Eq.  (8a),  we  find 

7  17  *  (C»!  a(°>)">S  ■  <«> 

If  the  frequency  mismatch  satisfy  the  inequality  Eq.  (22),  the  linear 
gain  expression  is  valid  for  z/2tt«  ^  Lb/Jt0. 
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For  the  low  gain  regime,  the  independent  non-dimensional 
parameters  are  C02a(0),  Aw/cj0  and  D02  .  The  appropriate  dependent 
variables  are  a(0)G  and  n*  In  the  low  gain  regime,  the  efficiency  is 
approximately 

ri  =  ^(1+  ^)(c-a<0')(a(0)G  )■  <24> 
IV.  Numerical  Results 

We  will  apply  our  •theory  to  illustrate  a  few  concepts  on  space 
charge  effects  and  filling  factors  using  a  practical  example 
in  the  low  gain  regime.  We  will  also  present  some  widely  applicable 
graphs  for  growth  rate  and  efficiency  in  the  high  gain  regime. 

a)  Low  Gain  Example 

The  FEL  experiment  we  have  chosen  consists  of  a  magnetic  wiggler 
with  a  wavelength  £0=  1  cm  (k0=  2tt  cm-1),  amplitude  B0=  325  gauss 
(A0  =  51.7  stat  volts),  an  interaction  length  L  =  3  m,  an 
electron  beam  with  energy  of  3  Mev  (y0=  7)  and  current  density  of 
16 A/cm2  (n0=  3.3  x  109em~3,  wb=  3.3  x  10”9sec_1).  The  radiation  is  at 
102  ym^36^.  Taking  filling  factors  to  be  unity,  we  obtain  £  =  6.5  *  10’3  and 

B0jl  =  4.33  x  10~3.  For  the  initial  example,  A(0)  is  chosen  to  achieve 

saturation  at  z=L  as  predicted  by  linear  theory^24).  For  the  parameters 
here,  A(0)  =  -  3.5  x  10"2stat  volts.  The  normalized  parameters  are 

C02a(0)  =  2.5  x  10"6  and  D02  =  2.85  x  10‘7. 

Figure  2  consists  of  plots  of  linear  gain  with  space  charge 
effects  (solid  curve),  02  =  0.76,  and  linear  gain  without  space  charge 

effects,  (dashed  curve),  92  =  0,  as  a  function  of  the  frequency  mismatch 

Au)/w0  at  z/2tt  =  300  (z  =  3m).  We  see  that  space  charge  effects  reduce  the 
linear  gain  by  about  20%  at  the  maximum.  The  dotted  curve  is  the  gain 
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0.4 


L  =  3m 


<  0.2 


FREQUENCY  MISMATCH - (10“3) 

COq 

Fig.  2  —  Plots  of  gain  as  a  function  of  frequency  mismatch,  Au/u0  at  3m.  The 
solid  curve  is  linear  gain  with  space  charge  effects,  0p2  -  0.76.  The  dashed  curve 
is  linear  gain  without  space  charge  effects,  0p2  =  0,  and  the  dotted  curve  is  the 
gain  at  3m  calculated  from  the  non-linear  self-consistent  equations  (8a-d). 

at  3m  calculated  from  the  non-linear  equations  (8a-d),  and  the  peak 
value  of  the  gain  has  shifted  to  a  larger  frequency  mismatch. 

If  the  frequency  match  is  small  relative  to  the  trapping  potential, 
as  in  Eq.  (20),  saturation  effects  become  important  when  the  electrons  travel 

roughly  half  of  a  bounce  length.  A  bounce  wavelength  is  the  distance 

electrons  travel  in  the  axial  direction  while  executing  a  bounce  period 

in  the  ponderomoti ve  well.  For  this  example,  — |  =  6  x  10" 3  and 

Y°trap 

L^/2  =  3.1m.  Therefore,  linear  theory  is  not  applicable  at  3m.  Plots 
of  linear  gain  with  space  charge  effects  (solid  curve),  linear  gain  with¬ 
out  space  charge  effects  (dashed  curve),  and  gain  from  non-linear  self- 
consistent  calculation  (dotted  curve)  as  a  function  of  axial  distance 


•  ,-v.  . 


for  Aw/u)0  =  -1.5  x  10'3  are  shown  in  Fig.  3  .  This  value  of 
frequency  mismatch  gives  the  maximum  linear  gain  at  z/2n  =  300  (z=3m) 
see  Fig.  2  .  The  linear  regime  of  the  non-linear  self-consistent 
formulation  is  z_<2m  for  this  example.  Finally,  we  will  show  that,  in 
the  linear  regime  of  the  non-linear  calculation,  the  gain  from  the  non¬ 
linear  calculation  is  in  good  agreement  with  the  linear  gain  expression 
as  shown  in  Fig.  4  calculated  at  z/2tt  =  150  (z  =  1.5m). 

To  summarize,  the  disagreement  between  linear  gain  and  gain  from 
non-linear  calculation  at  z  =  3m  is  that  the  initial  guess  of  CQ2  a ( 0 )  = 

2.5  x  10' 6  (A(0)  =  Al)  was  slightly  too  large.  Saturation  effects  have 
become  important  at  the  end  of  the  interaction  region.  As  CQ2  a(0)  or 
A(0)  decreases,  the  gain  from  non-linear  calculation  will  converge  to 
the  linear  gain  expression  as  shown  in  Fig.  5  .  Figure  5  contains 
curves  of  the  gain  from  the  non-linear  formulation  (solid  curves)  for  various 
values  of  C02a(0),  linear  gain  with  space  charge  effects  (dashed  curve), 
and  the  gain  from  non-linear  calculation  without  the  space  charge  potential 
and  C02a(0)  =  1.25  x  10”6  (A(0)  =  %  A^Hcircles)  as  a  function  of  frequency 
mismatch  at  z/2tt  =  300  (z  =  3m).  In  the  linear  regime  of  the  non-linear 
calculation,  the  space  charge  potential  decreases  the  gain,  agreeing 
with  the  linear  theory.  Figure  6  are  plots  of  efficiency  versus  frequency 
mismatch  at  z/2tt  =  300  (z  =  3m)  for  the  same  range  of  C0'a(0). 

The  effect  of  space  charge  on  gain  for  the  same  FEL  example  will 

be  summarized  below.  According  to  linear  theory  ,  the  space  charge  effect 

2 

is  only  moderately  important  for  9^  =  0.76.  We  found  that  the  effects  of 
space  charge  are  dependent  on  the  radiation  vector  potential  at  input, 

A(0).  Taking  A(0)  =  A^,  we  obtain  C02a(0)  =  2.5  x  10  and  g-z^jy  =  0.67. 

The  results  of  self-consistent  non-linear  calculations  found  the  space 
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. A: 


AXIAL  DISTANCE  z  (meters) 

Fig.  3  —  Plots  of  gain  as  a  function  of  axial  position  at  the  frequency  mismatch, 
Aco/oo0  *-1.5  X  10“3,  which  gives  the  largest  value  of  linear  gain.  The  solid  curve 
is  linear  gain  with  space  charge  effects,  6p2  =  0.76.  The  dashed  curve  is  linear 
gain  without  space  charge  effects,  d2  =  0,  and  the  dotted  curve  is  the  gain  from 
the  nonlinear  self  •consistent  equations  (8a-d).  Saturation  effects  are  important 


Fig.  4  —  Plots  of  gain  as  a  function  of  frequency  mismatch,  Acj/cj0,  at 
z  =  1.5m,  which  is  within  the  linear  regime  of  the  self-consistent  calcu¬ 
lation.  The  solid  curve  is  linear  gain  with  space  charge  effects,  6p2  = 
0.76.  The  dashed  curve  is  linear  gain  without  space  charge  effects,  6  2 
=  0.  The  dotted  curve  is  the  gain  from  non-linear  self-consistent  calcula¬ 
tions  (8a-d). 


Z/27T  =  300 
el  -  0.76 
D2  =  2.85x10 


C  2a{0)  =  6.25x10“7 


\V_  Co2a(0)  =  1.25xl0“6 

\Y\^  Co2a(0)  =  2.5X106 
^ Co2a(0)  =  5.0x  10® 


Fig.  5  —  Plots  of  gain  as  a  function  of  frequency  mismatch  at  z/2ir  =  300  (z=3m). 
The  solid  curve  is  gain  from  non-linear  self-consistent  formulation  for  various  val¬ 
ues  of  Co2a(0).  The  dashed  curve  is  linear  gain  with  space  charge  effects,  and  cir¬ 
cles  (o)  are  gain  from  non-linear  self-consistent  calculation  without  space  charge 
effects  with  C2a(0)  =  1.25  X  10'6. 


charge  effects  to  be  negligible.  Space  charge  effects  become  increasingly 
more  important  as  we  decrease  C02a(O)  or  A(0).  In  Fig.  5  ,  we  see  that 
the  gain  in  the  linear  regime  of  the  non-linear  calculation  without  the 
space  charge  effects  (circles  (o)  is  about  10%  larger  for  0^(0)  =  1.25x10" 

o  n  2 

(A(0)  =  1/2  A^)  and  _ =  1.3.  Figure  7  contains  curves  of  the  spatial 

evolution  of  the  amplitude  of  the  space  charge  potential  (dashed  curve)  and 
the  amplitude  of  the  ponderomotive  potential  (solid  curve)  for  Aoj/oj0  = 

-2.0  x  10'3,  C02a(0)  =  1.25  x  10'6(A(0)  =  1/2  Aj  and  002=  2.85  x  10'7. 

Notice  that  at  3m,  the  space  charge  potential  is  about  2/3  of  the  ponderomotive 
potential.  Thus  we  expect  the  space  charge  potential  to  have  a  moderate 
effect  on  the  gain. 

Figure  8  contains  plots  of  gain  with  space  charge  effects  (solid 
curve)  and  gain  without  space  charge  effects  (dashed  curve)  as  a  function 
of  axial  distance  for  =-2.0  x  10"!  C  2a(0)  =  1.25  x  lO’^AfO)  =  1/2  A.  ) 

0J0  0  L 

and  D02  =  2.85  x  10"?.  As  noted  in  Fig.  5,  the  space  charge  effects  decrease 
the  gain  in  the  linear  regime  of  the  self-consistent  calculation.  The 
collective  effects  of  the  space  charge  field,  however,  increase  the  gain 
in  the  non-linear  regime. 


Let  us  now  turn  our  attention  toward  the  effects  of  filling 
factors.  When  the  electromagnetic  filling  factor  Fgm  is  no  longer 
equal  to  unity,  the  appropriate  parameter  for  comparison  purposes  is 
efficiency  as  expressed  by  Eq.  (11a).  Plots  of  efficiency  versus 
axial  distance  for  various  combination  of  filling  factors  are  shown  in 
Fig.  9  for  =  -2.0  x  10~\  C02a(0)  =  1.25  x  10"6and  D02/Fsc=  2.85x  10"7. 
In  the  linear  regime,  z/2tt<300  the  efficiencies  calculated  with  different 
filling  factors  are  approximately  the  same.  In  the  non-linear  regime 
z/2tt>300,  efficiency  became  smaller  as  Fgm  decreases. 


AMPLITUDE  OF  POTENTIAL  (10 
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AXIAL  DISTANCE  zl2n  (102) 

Fig.  7  —  Plots  of  the  amplitude  of  the  space  charge  potential  (dashed  curve)  and  the 
amplitude  of  the  ponderomotive  potential  (solid  curve)  as  a  function  of  axial  position 
for  Aoj/cj0  =  -  2.0  X  10'3,Co2a(0)  =  1.25  X  10'6  (A(0)  =  1/2  AL)  and  DQ2  =  2.85 


(a(z)  -  a(0))/a(0) 


EFFICIENCY  (10’ 


=  -2.0X10"3 


Co2a(0)  =  1.25  x  10' 


D02  =  2.85x10 


^em  0.5, 
Fsc  =  0.95 


^em  Fsc  ^  ^ 

/  /  / 

/  /  / 

/  y 

// 

/■  /^Fem  =  0.1.  Fsc  =  0.95 

/// 

// 


AXIAL  DISTANCE  zilTT  (102) 

Fig.  9  —  Plots  of  efficiency  versus  axial  distance  for  various  combination  of  filling  factors 
with  A <0/wo  =  -  2.0  X  10'3,  C_2a(0)  =  1.25  X  10"6,  and  D^/F,,.  =  2.85  X  10’7 


In  the  low  gain  regime  the  non-dimensional  value  a(0)G  at  any 
axial  position  z  is  only  a  function  of  the  independent  parameters 
Cfl2a(0),  D02  and  Au)/co0.  A  set  of  three  curves  of  a(0)G  at  saturation  for 
different  non-dimensional  values  of  CQ2a(0)  is  shown  as  a  function  of 
frequency  mismatch,  Ad/u>0.  in  Fig.  10  for  D02  =  2.85  x  10“?. 
b)  High  Gain  Regime 

For  the  high  gain  regime,  we  illustrate  examples  where  collective 
effects  are  negligible.  Figure  11  is  a  plot  of  growth  rate  in  the  linear 
regime  of  the  non-linear  self-consistent  formulation  as  a  function  of 
frequency  mismatch  with  O02  =  5  x  10’7  and  various  values  of  CQ2.  The 
results  are  in  good  agreement  with  the  normalized  linear  growth  rate 
obtained  from  linear  dispersion  relation,  Eq.  (10).  The  saturation  ampli¬ 
tude  and  efficiency  as  a  function  of  frequency  mismatch  are  plotted  in 
Figures  12  and  13.  We  notice  that  the  magnitude  of  efficiency  as  well 
as  the  bandwidth  increases  as  C02 increases. 


GOq-O) 

CDq 


Fig.  10  —  Plots  of  a  (0)G  at  the  saturation  in  the  low  gain  regime  of  the  non-linear  t. 

self-consistent  calculations,  when  space  charge  effect  is  moderately  important,  D02 
=  2.85  X  10_7,asa  function  of  the  non-dimensional  parameters Co2a(0) and  Au/u 
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Fig.  12  —  Plots  of  efficiency  in  the  high  gain  regime  as  a  function  of  frequency  mismatch 
Aw/oj0  and  Ce2  for  a  case  when  the  collective  effect  is  negligible,  Dc2  *  5  X  10-7 
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Fig.  13  —  Plots  of  saturation  amplitude  of  the  radiation  vector  potential  in  the  high 
gain  regime  as  a  function  of  frequency  mismatch  Acj/o>0  and  Cc2  when  the  collec¬ 
tive  effect  is  negligible,  D02  *  5  X  10-7 
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